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Abstract 

The existence of higher derivative discontinuous solutions to a first order ordinary 
differential equation is shown to reveal a nonlinear SL(2,R) structure of analysis in the 
sense that a real variable t can now accomplish changes not only by linear translations 
t ^ t + h but also by inversions t 1/t. We show that the real number set has the 
structure of a positive Lebesgue measure Cantor set. We also present an extension of the 
Picard's theorem in this new light. 
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1 Introduction 



The most basic ingredient of real analysis is of course a variable, living in the real number line, 
which is assumed to undergo changes by linear translations. The framework of ordinary analysis 
is therefore essentially linear, the relevant concepts of limit, differentiability etc being formally 
defined to satisfy the simple transformation properties of the linear translation group. Notice, 
for example, that the limit t — )■ to means the linear, continuous decrease ( as an action of the 
translation group) of the distance |t— to| to zero and so on. The application of this linear analysis 
in more complex dynamical problems, for instance, necessitates explicit breaking of this linear 
structure by incorporating nonlinear interactions (terms) in the underlying differential (and/or 
integral) equations. The recent investigations, however, reveals the surprising existence of some 
nonlinear dynamical structures right at the heart of the ordinary analysis P, [s], 0, Is], thereby 
extending the framework of the linear ordinary analysis to the scale free SL(2,R) analysis. This 
extension is based on the realization that the simplest linear, scale free, differential equation, 
viz., 

4=...(i)=i (1) 



admits a generalized class of nonlinear solutions o: 



the form Tg = Ts{l+0{ri'^")), having discon- 



tinuities in the 2"th order derivatives at t = 1 jsj. Here, = t is the standard solution and 
t ^ t/to = 1 + f] denotes a rescaled variable close to 1. The discontinuity in the 2"th order 
derivative stems from the rescaling invariance of eq(Il]) which allows one to inject an infinite 
sequence of arbitrary, infinitesimal scale factors into the standard solution Tg. Consequently, 
the small scale variable t] is raised to a random infinitesimal variable, residing in a nonstandard 
real line R rather than in the ordinary real line R. Indeed, let t± = 1 ± rj, rj > 0. Then 
obviously Tg- = Ts(t_) = 1 —rj. Now, it is easy to verify that this exact solution can as well be 
derived starting from an initial approximation and then obtaining self similar (multiplicative) 
correction factors recursively, viz., r, = (l/t+)ri when the correction factor satisfies the self 
similar equation 

'-It-- 

in the smaller scale variable ti± = 1 ± r^^ in the neighbourhood of 1, and so on, since t_ = 
(i+r]){i+if ){i+ri'i) ■ '^^^ freedom of a residual ^ rescaling in eq(|2]) of the form ao^i- = 't'i_ = 



^The said rescaling is residual in the following sense. The trivial rescaling t t/to shifts a variable near to 



1 — rj', r]' = ao{r]'^ — eo/ao), ao = 1 + eo, Cq being an arbitrarily small real number drawn 
at random from the interval rj^ < to << f], however, induces a small scale modulations ( 
fluctuations) in t'i_ and so on to the higher order iterates generated self similarly from 
eq([2]) over smaller and smaller scales. The new, nontrivial solution jsj] is thus obtained as 

^^=^117^' ^=nc(o) (3) 

Clearly, the smaller scales reveals a nonlinear structure in the neighbourhood of a real number. 
Notice, for example, that r]'^ = an{o!^n^\fi^^^ — enjo-n)-, representing higher order uncertain- 
ties and antn+ 7^ t'^+- The second inequality tells the nontriviality of the residual rescaling 
symmetry, which is responsible for the higher derivative discontinuity of the new solution Tg at 
t = 1, provided, of course, the rescaling is performed an inflnite number of times ( for a proof 
see Appendix). Notice also that Tg is 2nd derivative discontinuous if the nontrivial rescalings 
are injected at the 2nd and an inflnite sequence of subsequent iterates, 4th derivative discon- 
tinuous when the freedom of this rescaling is used for the flrst time at the 3rd, instead of 2nd, 
iterate and so on. In general, the new solution would have discontinuous 2"^^ derivative if the 
nontrivial rescaling is utilized at the nth iterates (the standard solution is recovered when 
this rescaling is postponed till the inflnite number of iterations 5| (c.f. Appendix)). 

In view of this class of discontinuous solutions to the simplest scale invariant linear equation 
([T]), the ordinary real number system R needs to be extended, accommodating this small scale 
nondifferentiable structure, to a nonstandard -like real number system R. Notice that the 
generalized solution of eq([T]) can indeed be written as Tg = t{l ±e(j){ti)), (pit) = Tiy(t)/t, ti = et 
where = 0, r/v being the new solution ([3]). Consequently, every ordinary real number, for 

instance, to gets extended to a set of the form (to)e = {t '■ to—e<t< to+e, e 7^ 0, O(e^) = 0, t ^ 
to} CR where e denotes the irreducible uncertainty even when the real number is evaluated 



or 



with an "inflnitely precise measurement". Clearly, the inclusion . . . (to)e2 C (to)e is satisfled 
higher precision evaluation of to- We note that the generalized solution, as it is obtained in 
is constructed purely as a function deflned in R and hence the extended set R constructed out 
of the generalized solution coincides exactly with R viz., R= R. We remark that this feature 
of our extension distinguishes it from the work of Robinson Q]. In the traditional approach 
of nonstandard analysis, inflnitesimals are deflned as an equivalence class of sequences of real 
numbers. As a result, inflnitesimals, being constructed out of real numbers, can be thought of 



to one near 1. The present rescaling reveals, the so far hidden, freedom of injecting a small scale uncertainty in 
the neighbourhood of a point. 



as extraneous to the set of ordinary real number set. Our results, on the other hand, suggests 
that infinitesimals are indeed members of the real number set R and so might instead be defined 



formally even in the real number set R. This approach is being developed recently by Bose 
In view of this existence of infinitesimals in R, with random properties, the new solution G R 
can as well be interpreted to reveal a new mode of change that is available to a (real) variable 
t, viz., t can change near an infinitesimal neighbourhood of t = 1 by (random) inversions: 
t_ —7- tl^ = t+, t± = 1 ± ri, for an infinitesimal r], defined by 7^ but 0(?7^) = 0, besides the 
ordinary translations over the ordinary real number line. Notice that the scale free equation ([T]) 
is invariant under the transformation t— T-l/t, r— 7>l/r and so t_ — )■ tZ^ = t+ indeed stands 
for the generalized solutions considered here. The possibility of (random) inversions breaks the 
exact determinacy of a real variable. Further, the scale invariance tells that the small scale 
fiuctuating behaviour of a real variable ( over rj) is reproduced self similarly over the smaller 
nonlinear scales {rj'^"), though the self similarity is respected in a statistical sense. Accordingly, 
the linear framework of the ordinary analysis gets extended to a SL(2,R) analysis. One therefore 
infers that the real number system R, as it is ordinarily understood in the context of the linear 
(ordinary) analysis, enjoys a host of richer "dynamical" properties which get revealed in the 
SL(2,R) analysis. We have already discussed some of these new features of this analysis in 4|. 
Here, we present an analysis of the Picard's theorem 8|] in the light of the new generalized 
solutions of eq(IT]). To this end we first discuss how the arguments used in the simple equation 
([I]) can be extended to more general equations (Sec. 2). Next we discuss some salient features of 
the computational model introduced in 5j and also show that R has the structure of a positive 
Lebesgue measure Cantor set The (box counting) dimension of a point in R ( hence in 
R ) has the value z/, where u = ^ is the golden mean (Sec. 3). Before stating the Picard's 
theorem in this new light, we present an extension of the ordinary Riemann integration on R 
(Sec.4). 

2 General case 

To prepare for a discussion of the Picard's theorem, let us clarify further the set up of the scale 
free SL(2,R) analysis. To this end, we consider the more general linear ODE of the form 

/(*..). .(1)^.0 (4) 
Let f{t) be the corresponding standard solution. This solution is exactly determinable, in 



principle, in the linear analysis. In the extended framework this, however, corresponds to the 
zeroth order solution only. Based on this solution, one can, however, generate a more complex 
solution in the following way. Following the steps outlined in Sec.l, let us write the generalized 
solution close to t = 1 (more precisely in a Ihs neighbourhood of 1) as Tg-{tg_) = r(t_) = 
f(l/t+)ri. Then it is easy to verify that t'__ would satisfy an equation of the form 



on the smaller nonlinear scale ln(l — r/^) where /'(ti„, rl) = (*+-^(* ) ^ ^|iA+.'r(iA+))) ^ Clearly, 



the generalized solution Tg would belong to the higher order nondifferentiable classes of solutions 
when the infinite sequence of nontrivially rescaled (equivalently, random) variables t'^_ are 
incorporated as indicated above. 

Notice that the exact self similarity on the smaller (nonrandom) scales ln(l — rj'^) and 
etc is obtained only for the simplest equation ([T]). For, a t— dependence either of the form 
f{t) or f{t,x) breaks this exact self similarity. Nevertheless, an approximate self similarity is 
maintained on scales randomised by the nonlinear t-dependence, besides the above mentioned 
freedom of rescaling, when the random scale is defined, for instance, by t[_ = ^ 1 — 
fif]'^ where n = f'{ti_,T'_)/f{ti_,T'_) ^ 1, which remains almost constant over the scale ?7^. 
Note that, randomisation of scales by explicit nonlinearity in the ODE already simulates the 
residual rescaling freedom of eq([l]) . The emergence of random behaviour in dynamical systems 
(for instance, one dimensional maps, higher order nonlinear ODEs etc) from truncation errors 
are well known in deterministic chaos. Here, we point out equivalent behaviour even in one 
dimensional linear ODEs as an effect of infinitesimal nonlinear scales (elements) in real number 
system. As examples we consider two simple equations in the following. 

Case 1. Modulated exponential 



Small scale modulations in the ordinary exponential function defined by 

dt ^ ^ 

is expected in the present formalism. One verifies that the generalized solution near t = 1 is 
obtained as 




(5) 



2 

where ti+ = 1 + — t-) = 1 + so that /i = 1/(1 — r]), and the equahty in ([7]) is up 

to a multiphcative factor. In a computational problem (truncated) nonlinear terms in /i would 
build up over time to introduce small scale random oscillations in the ordinary exponential. 
Such behaviour would of course be common to any deterministic functions of linear analysis. 



Case 2. Quadratic nonlinearity 
Let us consider the equation 

t% = r\ r(l) = 1 (8) 

having the exact solution tq = 1/(1 — Int). Writing Tg{tg^) = ro(t_) = ro(l/t+)ri_, as usual, 
we get 

dlnri_ 
dh^ = 

where /i = i-int' ~ point out that an equation of the form ([9]) is truly 

nonlinear in the sense that nonlinear random scale t[_ = ^ 1 — ^rf is itself determined 
by the "unknown" function ri_ through /x and hence does not lead to an exact solution in the 
ordinary sense. 

To conclude this section, we remark that the infinitesimal fluctuations induced to the solu- 
tions to eq(jll) because of SL(2,R) inversions would have a universal feature determined by the 
functional dependence f(t,T) close to t = 1. 

3 Infinitesimals, computation and Cantor sets 

In ref 5| we presented a computational model of infinitesimally small real numbers. Here we 
explain the relationship of inversions and infinitesimally small numbers in the context of this 
model. Recall that every ordinary real number t is replaced by = t(l ± e0(ti)) so that 
essentially corresponds to a continuum set (t)^. Now in a computational problem, a real 
number is only treated as a finite precision (decimal/binary) representation. Consequently, the 
real number set is covered by a countable collection of disjoint open intervals of the form t^, 



where e now denotes the finite precision with which the real numbers are evaluated (calculated). 
Notice that e can have values, in decimal representation, of the form 0.5 (for integers, so that 1, 
for instance, stands actually for the open interval (0.5,1.5) and so on), 0.05, 0.005, 0.0005,. . .. As 
a result, in any computation the real number set, although conceived as a connected, continuum 
set (without any gap), is actually realized as a totally disconnected, countable set (equivalent 
to the set of rationals). As one improves upon the accuracy and approaches to the infinite 
precision the connected continuum set of linear analysis is thought to have been recovered only 
in the limit. In view of the scale free extension, one, however, meets with an obstacle in the form 
of irreducible infinitesimal uncertainties, thereby realizing the continuum structure of the real 
number set but for the total disconnectedness. In fact, the extended real number set acquires 
the structure of a Cantor set with a positive Lebesgue measure Gj. 

Lemma: The set R has the box counting dimension a, a = 1 + ln(l + e(f))/ Int. Further, a 
has the asymptotic limits limf^o+ cr = z/ and limj^oo c = 1 + u where u = . 

Proof: In view of the structure of new solution presented in Sec.l, R is an extension of R 
in the following sense. Let A^^o^ be the number of open balls (intervals) of radius 6 covering an 
open interval of R. The box counting dimension of R is lim5_>o+ which, in fact, equals 
1. To cover R the set {Nqs}, however, is not sufficient. One, in fact, needs to subdivide each 
of subintervals of size 6 into smaller subdivisions of size 5^ and this process of subdivisions 
is continued recursively over smaller and smaller scales 6"^, 6^ . . .. In each step R is covered 
better but for a tiny residual part, which could only be measured (detected) looking at the 
next level of finer scale. Accordingly, these subdivisions over finer scales facilitates one to 
measure the required "length" of R more and more accurately only in a progressive manner^ 
Let, in the first step of subdivisions, the fractional increase in the number of subdivisions be 
Nis2/Nqs{= a < 1), where Nfg is the number of subdivisions of scale 5^ needed to cover the 
boundary "points" of the original interval. Notice that the interior points in the open interval 
concerned are already covered well by the zeroth level covering balls and so finer scale covering 
balls are irrelevant for these points. Continuing this process sequentially over smaller scales 
(as indicated above) one then gets, for instance, as the fractional increase N2si/Nip and so on. 
Because of the (relative) self similarity of the subdivisions one concludes N2S'i/Nis2 = Nis2/Nqs 
etc. Consequently, total number of subdivisions of different sizes needed to cover the said 
interval of R is A^ = ^ N^s^" = Nqs{1 + A + + . . .). Hence the box counting dimension of 



^The concept of a perfect (or total) covering (and/ or infinite precision) is, however, illusory. What is more 
meaningful is the continual progress towards perfection 0, Q . 



R, evaluated relative to scale S, is obtained as a = lim5_>.o+ = 1 + 

Let t G R. Then the "length" of the open interval (0, t) would be given by where t & R. 
Again, in view of of the new solution we write t = t(l + ecj)). Comparing the two we deduce 
the first result. The value of A is thus determined by the nontrivial factor in t. Its exact value 
is, however, unimportant for our purpose. That it is nonzero is actually more significant. 

To prove the remaining part of the lemma, we note that \mit^Q+ = linit-i^oo(l~e( i^7Pt)t(^i) + 
higher order terms) = limj-i^oo(l — ^( in^-i + + • • •) denotes the new 

solution of ([3])). As t^^ approaches to oo through higher and higher order scales, there exists a 
sufficiently large interval in i?, when 7r(t^^) = j^-pr ~ ~ e^^. Accordingly, in this interval 
limt_,o+ C7 = limt-i^oo(l/(l + r(ti))) = \imt-i^^{l/ {1 + (1/(1 + e\^)T{t2))) where = eH 
and so on. Notice that the r.h.s of the first equality mimics exactly an application of inversion 
on the first two terms of the logarithmic expansion of a. The desired limit now follows from 
the continued fraction expansion of the golden mean v. Similar arguments also hold for the 
other limit when one makes use instead 0(t) = tr(t^^), r(t^^) = l/r(ti) [3]. □ 

Remark 1 To see the origin of inversion, let us reconsider the computational model in- 
troduced above. Let a computation could only distinguish integer numbers, so that 1 = (1).5, 
for instance. Consequently, (1)^ = 1 + e/x, = /i_ U /i_ = ( — 1,0), /i+ = (0,1), e = 0.5. 
The number 1, therefore, corresponds to any number drawn perfectly at random from (0.5,1.5) 
(the assumption of randomness is codified in the residual rescaling symmetry c.f. Sec. 1). 
This is realized in the above by letting e instead to lie in the interval (—0.5,0) U (0,0.5) and 
/i = ~0(1). The irreducible fiuctuations as encoded in the new solution now tells that a 
point in (-0.5,0) could fiuctuate to a point in (0,0.5) (by inversion, viz., by scaling and fiipping 
of sign) and vice versa, instead of pure translations, as assumed in ordinary analysis. Further, 
the especial role of inversions in a Cantor set is also intuitively clear because a point in a Cantor 
set can change to another only by an inversion (i.e., by a discrete jump) because of the gaps 
(voids) separating the Cantor points 4|. 

Remark 2 To visualize the Cantor set like structure of (0)^, we note that (0)^ = (— e, Oq) U 
(0)^2 U (Oq, e), where Oq denotes the zero (0) at level 0, and so on recursively over smaller scales. 
Accordingly, in this extended framework, the number 0, for instance, is represented by the set 
U[f {(— e^", 0„)U(0)^2"+i U(0„, e^")}, 0„ being the nth level zero. Because of the natural inclusion 



On D 0„+i, the set {0„} represents a finer and finer realization of the originally coarse grained 
zero (0) as the accuracy of the computation is increased. In linear analysis, however, these fine 
structures are ignored so that the ordinary singleton {0} is reproduced. Notice that the usual 
singletons {to} in the reduced (coarse grained) set R are reinserted once the nontrivial SL(2,R) 
generator is frozen. 

Remark 3 The significance of vr(t~^) in the prime number theorem and other number 
theoretic results will be considered elsewhere. 

4 The Picard's theorem 

In view of the Cantor set like structure of R (every point of R is actually got replaced by a 
Cantor set, c.f., remark 2), the definition of the ordinary (Riemann) integration needs to be 
extended. Notice that the Cantor set of (1)^, for instance, is written as (1)^ = U^[(l„ — e^", 1„ + 
e'") - {In}] (c.f., remark 2) of length 2 - (e'" - e'"^')] = 2 e^" which is essentially 

zero, being of higher order infinitesimal, relative to the first order infinitesimal e. Notice that 
the length of voids in (1)^ is 2 J^i^"^" " e^"^^)) = 2e, which equals to the length of (1 — e, 1 + e) 
in R, so that the Lebesgue measure of this Cantor set is effectively zero. 

Notice that our iteration process, leading to the new solution of eq(IT]), not only reveals a 
multiplicative structure of the nontrivial neighbourhood of 1, but also factors the differential 
operator in ([T]) into a countably infinite set of self similar operators over nonlinear scales t'^. Ac- 
cordingly, the corresponding integration measure defined on the said Cantor set, incorporating 
random inversions, is defined by the following replacement 

/ dlnt.^^R/ dlnt.^V / dlnt^^ (10) 

Jl-rj Jl-n ^ Jl 

where rj is an infinitesimal variable, so that the new solution of eq([I]) Tg = CYl jr—, C = 

n 

n^n+(0) pI is retrieved by direct integration from eq([T]). Notice that 1 in the integral of the 
r.h.s. of eqf llOp corresponds actually to 1„ with the definition that In 1„ = for each n. More 
generally, the extended integral for a function / in R is defined by 

ff{t)dt= l"fXt)dt + eJ2 r [f{t'-') + Tgt^]d\nt'n (11) 

J Sl J Ob J \ 

where we neglect higher order infinitesimals, and use rescaling invariance t ^tjaoi the loga- 
rithmic differential. We also make use of Tg = t^, dr^ = Tg'^dlnt'^ and t = t + erg. Notice 



that the second term in the bracketed integral arises from the first order Taylor's expansion of 
/(t) and the extended integral can be considered as the anti-derivative of the ODE 



f = /(t) (12) 

in R (c.f. eq(jlj)). As an example, the modulated exponential ([7]) is recovered when eq(l6]) is 
integrated following the generalized integral ( ITTj) with f{t) = 1. 

Remark 4 Clearly, the infinite sum in ( ITT]) is uniformly convergent in every closed interval 
near t = 1 for continuous df/dt. Let us also point out that the infinite set of rescaled nonlinear 
variables t'„ are intrinsically randomised, as indicated in relation to eq(|l]), implicitly by f{t) 
itself. Consequently, the extension of integral as defined in eq f|T2|) should be interpreted as a 
generalized integral over a set with small scale (infinitesimal) random elements. We make a 
more rigorous treatment of this generalized integral elsewhere. 

Proceeding at this heuristic level (i.e., forgetting the randomness so that all the variables 
t and are well defined functions of the real variable rj, < 1] << 1) we now reexamine the 
Picard's theorem in this new light. 

Let us consider eq(jl]), viz., 

^ = /(t,r),r(l)=ro (13) 

where / is assumed to be a C°° function, for simplicity. The Picard's theorem then guarantees 
a unique solution in the neighbourhood of t = 1. It also tells that the corresponding solution 
is C°°. The proof requires one to convert the I VP f|T3l) to the equivalent integral equation 

Inr = lnro + ^ /(t,r)dlnt (14) 
and then to construct a sequence of approximations r„(t) satisfying 

lnr„ = Inro + ^ /(t,r„„i)dlnt (15) 

converging uniformly in a closed interval within the said neighbourhood to the required solution. 
One verifies that the proof of the Picard's theorem applies also to the generalized integral (fTT]l 
without major modifications. Consequently, we state the modified Picard's theorem in the 
SL(2,R) analysis as follows 



Theorem Let /(t,r) be C^""! in R. Then the ODE 



^ = /(t,r),r(l) = ro (16) 
has a unique solution in a suitable neighbourhood of t = 1. 

Remark 5 The function / in the above theorem inherits the differentiable structure of t. 
If instead / G C™, m < 2" — 1 then the solution would obviously be only C™. 

Appendix 

The new solution tn reduces to the standard solution if the residual rescalings are performed 



only up to a finite number of iterations, 
so that t'(„+i)+ = 1 + r]n+i, Vn+i = a^r]'^ 
then lead to a factor 1 — vfn+i) '^n, viz. 



et the rescaling be terminated at the nth iteration, 
sj]. Subsequent iterations (without rescalings) would 



^N- = C^-j- . . . -^(1 - V(n+1)) 
''+ ''2+ ''n+ 

Consequently, the uncertainty (or twist, to put it in a pictorial way ) introduced in the solution 
through each nontrivial rescaling and subsequent inversion would unwind gradually, so to speak, 
leading to the ordinary (nonrandom) Tg (multiplied by n irrelevant scale factors). Things, 
however, change drastically if rescalings are performed an infinite number of times, instead. 
Indeed, the statement n tends to infinity means, of course, that n is becoming larger and 
larger, signifying an unending process. As a result, the unwinding can never be initiated in the 
iteration process for an infinitely large n. 
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